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1. Introduction and Summary 

The acceleration of the universe was discovered 10 years ago and its explanation is one 
of the most secret enigmas of current theoretical physics. According to and the other 
experiments that map the anisotropics of the cosmic microwave background, the approx- 
imately 76% of the energy content of the universe is not dark or luminous matter but it 
is instead a mysterious form of dark energy that is exotic, invisible and unclustered. In 
order to explain the origin of this form of matter three main classes of models for this 
acceleration exist: 

1. A cosmological constant A 

2. Dark Energy 

3. Modified Gravity 

The cosmological constant is the most obvious explanation of the acceleration of the uni- 
verse |U . However A suffers from well known cosmological problem that requires extreme 
fine-tuning. 

The second class of models postulates the existence of a dark energy fluid with equation 
of state P ~ — p where p and P are the energy density and pressure of the fluid 

The last class of models known as Extended Theories of Gravity (ETG) corresponds 
to the modification of the action of the gravitational fields. These theories are based on 
idea of extension of the Einstein Hilbert action by adding higher order curvature invariants 
such as R 2 , R^ v R^ v .... and/or minimally or non- minimally coupled scalar fields to the 
dynamics . In this paper we consider the models corresponding to f(R) theories that are 
general functions of the Riemann scalar R. 

As we said above f{R) theories should mainly explain the mysterious acceleration of 
the Universe and hence provide modification of the Einstein theory on large scales. On 
the other hand they do not solve the second fundamental problem considering theory of 

1 For review and extensive list of references, see |(| [?], ^, ^| pxft . 
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gravity which is the fact that the quantum gravity is non-renormalizable theory. In order 
to improve this situation new formulation of the quantum gravity theory was proposed 
recently in series of papers 11, [l^, |||, [l3|] 2 . The key idea of this proposal is to abandon 
the fundamental role of the local Lorentz invariance and assume instead that this appears 
only at low energies as an approximative symmetry. The breaking of Lorentz invariance is 
achieved by introducing additional structure to the space-time that is a preferred foliation 
by 3— dimensional space-like surfaces and then splitting of the coordinates into space and 
time. This procedure allows to complete the Einstein's general relativity with higher spatial 
derivatives of the metric and then improve the UV behavior of the graviton propagator 
and makes the theory renormalizable by power-counting. 

Even if the Hofava-Lifshitz theory has many attractive properties investigations of 
various aspects of Hofava gravity have started to reveal some potentially troublesome 
features. The first one is considering the detailed balance condition that strongly restricts 
the form of the potential terms in the gravity action. However the detailed investigation of 
the Hofava-Lifshitz theory performed in past few months strongly suggests that this theory 
has much better behavior when we abandon the condition of the detailed balance 3 . 

The second fundamental ingredient of the original Hofava-Lifshitz was the projectability 
condition that says the lapse function in ADM formulation of gravity depends on time only. 
Even if for the most important vacuum and cosmological solutions Einstein equations 
can be put into the form of constant lapse function it seems that the consistency of the 
Hofava-Lifshitz theory requires the projectability condition. In fact, the investigation of 
the algebra of constraints of the Hofava-Lifshitz theory without projectability condition 
performed in paper p3| implies some pathological behavior of this theory. On the other 
hand it was already shown in |l2j that the algebra of constraints of Hofava-Lifshitz theory 
with projectability condition is closed. However it is important to stress that we should 
be very careful considering the definite conclusion about the consistency /non-consistency 
of Hofava-Lifshitz theory with/without projectability condition as was discussed recently 
in |l5|, [?2| in relation with an another issue of Hofava-Lifshitz theory that is spin-0 scalar 
graviton with addition to the standard spin-2 tensor graviton. The study of the properties 
of this mode showed that it persists down to low energies and exhibits the pathological 
behavior that should be disaster for the consistency of Hofava-Lifshitz theory. Explicitly, 



it was shown in 32] that this mode is strongly coupled at energies above a very low energy 
scale and it suffers from fast instabilities. It turns out that these properties appear both in 
the non-projectable and projectable cases. On the other hand it was argued in very recent 
paper [^] that it is possible to extend the original Hofava-Lifshitz action in such a way 
that these dangerous properties are absent. 

Despite of these open problems it is clear that Hofava-Lifshitz theory is very interesting 
and stimulating proposal that should be extended in many directions. Our goal is to 
generalize this theory in the similar way as the ordinary Einstein gravity is extended to 



2 For detailed study of this proposal, see 14 
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f(R) theories. We started this program in our previous paper [28] where we formulated the 
Hamiltonian form of f(R) Hofava-Lifshitz theory based on the condition of the detailed 
balance. We considered one explicit example of f{R) Hofava-Lifshitz theory with the 
function f(x) = y/x and we found the Lagrangian version of this theory. In this paper we 
generalize this analysis to the case of arbitrary f(R) Hamiltonian form of Hofava-Lifshitz 
theory. We show that when we extend the original Hamiltonian of f(R) Hofava-Lifshitz 
theory by introducing two non-propagating degrees of freedom we can find the Hamiltonian 
in such a form that allows us to perform the Legendre transformation from the Hamiltonian 
to the Lagrangian description. Then solving the equations of motion for these auxiliary 
fields we finally find a new form of f(R) theories that are invariant under foliation of 
preserving diffeomorphism [|i~2^ [. 

It is also possible to introduce the Hofava-Lifshitz theories of gravity in an alternative 
way when we replace 4-dimensional Ricci scalar R with combination R — ► KijQ^ 1 K^i + . . . 
known from the linear Hofava-Lifshitz theory of gravity. The resulting theory is invariant 
under foliation preserving diffeomorphism. 

We hope that this short note could be the starting point in the study of the new form 
of f(R) theories of gravity. It will be certainly very interesting to study the equations of 
motion that follow from these actions and find their possible cosmological solutions. It 
would be also interesting to perform the analysis of the fluctuation modes around classical 
background and study their physical properties. We hope to return to these problems in 
future. 

2. The First Version of f(R) Hofava-Lifshitz Theory of Gravity 

In this section we consider f(R) Hofava-Lifshitz theory of gravity that was introduced in 
p8]|. These theories are defined by Hamiltonian 



H = J d D x(N(t)H (t,x)+N i (t,x)H l (t,x)) , (2.1) 

where TL l is the generator of the spatial diffeomorphism 

H\t,x) = -2V i 7r ji (x) (2.2) 

that takes the same form as in standard Einstein theory of gravity. On the other hand the 
Hamiltonian constraint TCq of f(R) Hofava-Lifshitz TCq is defined as 

Wo(*,x) = K 2 ^f (Q^^g ijkl Q kl ^j , (2.3) 

where k is the coupling constant of the theory, g = det gij and where / is an arbitrary 
function. Further, the functions Q lJ and Q^ lJ are defined as 

Qij = i7T ij + ^g E ij ; q\h = _ i7r ij + ^ E ij ) (2.4) 

where g^, ,i,j = l,...,D are components of metric and tt 13 are conjugate momenta. These 
canonical variables have non-zero Poisson brackets 

{^(x),vr fc '(y)} = i( ( 5f4+^)5(x-y) . (2.5) 
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Finally we explain the meaning of the objects E 13 . In the original version of Hof ava-Lifshitz 
theory based on the the condition of the detailed balance they are given by variation of the 
D-dimensional functional W = f d D Ji.^/gw(g,dg) 

1 5W 

V^(x) = --^. (2.6) 

2 6gij{x) 

However we can clearly abandon this requirement and consider the generalized Hofava- 
Lifshitz theories where the Hamiltonian constraint 7i$ takes the form 

n = K 2 ^f (^7T ij g m 7r kl + E i3 G ljkl E kl + , (2.7) 

where . . . denote the contribution from additional terms that violet the detailed balance 
condition and that are crucial for reproducing General relativity at low energy regime. 
However for simplicity of the notation we consider the first form of the Hamiltonian ( |2.3| ) 
keeping in mind that the analysis below is valid for the general case as well. Finally the 
"metric on the space of metric", Q % i kl is defined as 

gijkl = ^ {gikg jl + gilgjk _ V ^H } (2 g) 

with A an arbitrary real constant. The inverse Qijki takes the form 

Gijkl = -^{9ik9jl + 9il9jk) — ^9ij9kl ; (2-9) 
where A = D ^_ l - Note that (|2~8|) together with (|2~9|) obey the relation 4 

QijmnQ mnhl = Usty + S^Sj) . (2.10) 



We should stress one important point. The Hamiltonian ( [2.1| ) contains the lapse function 
that depends on t only. This property is known as projectability condition. It is still an 
open problem whether this condition should be imposed or whether we should consider the 
lapse function with the full space and time dependence. In this section we consider the 
theories that obey projectability condition. When this condition is imposed the functions 
N(t) and Ni(t, x) have conjugate momenta 7T7v,7r J (x) that form the primary constraints of 
the theory 

tin ~ , vr^x) w . (2.11) 

Since these constraints have to be preserved during the time evolution of the system we 
find 

d t n N = {tt n , H} = J d D xH (x) = T T « . (2.12) 



4 We use the terminology introduced in and that we review there. In case of relativistic theory, the 
full diffeomorphism invariance fixes the value of A uniquely to equal A = 1. In this case the object Qijkl is 
known as the "De Witt metric". We use this terminology to more general case when A is not necessarily 
equal to one. 
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Note that the global constraint TV ~ is fundamentally different from the local constraint 
TCq (x) that is imposed in case when the theory does not respect the projectability 
condition. Further, the consistency of the constraints 7T*(x) ~ imply the secondary 
constraints 

a t vr'(x) = (tt^x), H) = -H*(x) « . (2.13) 
It is natural to introduce the smeared form of the spatial diffeomorphism constraint 



d"xCi(x)?f (x) . (2.14) 

The next goal is to calculate the Poisson bracket of constraints (2.12) and ( |2.14| ). These 
brackets have been determined in p8| with following result 



{T T ,T T } = , 
{T s (C),T r } = , 
{T S (C),T S (£)} = T s (Cd i e-Cd l C k ) . 

(2.15) 

In other words the constraint algebra of the gravity action that is invariant under foliation 
preserving diffeomorphism is closed. This is a crucial simplification with respect to theories 
that do not respect the projectability condition. More precisely, the algebra of constraints 
of the Hofava-Lifshitz theory without projectability condition was calculated in [43| with 
the result that suggests inconsistency of this theory. 

2.1 Lagrangian Formulation 

In this section we formulate the Lagrangian formalism for the f(R) theories defined by the 
Hamiltonian ([2.1|), In order to this we argue that the Hamiltonian constraint TLo defined 
in (12. II) can be written as 



W = K 2 ^/g 



B(o^-g ijM Q kl -A \ +f{A) 



+ v a Pa + v b Pb , 

(2.16) 

where we introduced two new fields A(x), B(x) with conjugate momenta Pa(x), -Pb(x) and 
corresponding Poisson brackets 

{A(x),P A (y)} = S(x-y) , {fl(x),P B (y)} = $(x-y) . (2.17) 

Looking on the form of the Hamiltonian density ( p. 16 ) we see that it is natural to interpret 
Pa (x) , Pb (x) as the primary constraints of the theory 

P A (x)«0, Pb(x)«0. (2.18) 

Then the requirement that the constraints ( |2.18| ) are preserved during the time evolution 
implies an existence of the secondary constraints Gvl(x), C?b(x) since 

8 t P A = {P A ,H} = k 2 ^(B - /'(A)) « k 2 ^7jG a « , 

d t P B = {Pb,H} = -K 2 yfg{Q^-G ijkl Q kl - A) = -k 2 ^G b ~ , 

(2.19) 
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where H is given in ( [2,1| ) with Hq defined in ( 2.16| ). Let us now calculate the Poisson 
brackets between P a ,Pb and G a ,Gb- It is a straightforward exercise to calculate them 
and we find 

{P A ( X ),G B (y)} = <5(x-y) , 
{P A (x),G A (y)} = f"(A)6(x-y) , 
{P B (x),G A (y)} = -<5(x-y) , 
{P B (x),G B (y)} = . 

(2.20) 

These results imply that P A , Pb,G a ,Gb form the second class constraints. As a con- 
sequence the requirement of the preservation of the constraints G a ,Gb under the time 
evolution fixes v a ,vb uniquely. Since these constraints are second class they strongly van- 
ish and can be explicitly solved with the result 

1 

-< 

9 

Then if we insert these results into the Hamiltonian density given in ( 2.16| ) we find that it 
takes the form 

/ . . . 1 . A 

(2.22) 



B = f'(A) , A = Q m -g m Q kl 



(2.21) 



n = K 2 ^gf (Q^G m Q kl ] 

that coincides with ( |2.1| ). As the final point note that for the system with the second class 
constraints we have to introduce the Dirac brackets instead of Poisson brackets 

{^(x),^(y)} D = {^(x),vr fc '(y)}- 

- J dx / dy'{^(x),cI> / (x / )}^ J (x / ,y / ){$j(y , ),vr fcZ (y)} , 

(2.23) 

where I = 1, 2, 3, 4 is a collection of the second class constraints <3?/ = (P A , G A , Pb, Gb) 
and where O is inverse of the matrix of Poisson brackets ( 2.20 ). Explicitly 



7J (x,y) ^{$ / (x),$ J (y)} 



/ /" 1 \ 
-/" 1 
-10 
\ -1 0/ 



<5(x 



(2.24) 



so that 



/J (x,y) 



<5(x 



(2.25) 



/o 1 \ 
0-1 
1 -/" 
V 1 /" / 

However due to the structure of the matrix Q IJ and the fact that g^j and ir 1 ^ have van- 
ishing Poisson brackets with P a ,Pb we find that the Dirac brackets coincide with Poisson 
brackets. 
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In summary, we claim that the Hamiltonian density that is useful for the developing 
of the Lagrangian formalism takes the form 



H = 2. Y,-V + Nk 2 ^ 



B[Qt ij ~g m Q Kl -A) +/(.!) 



1 



+ v a Pa + v b Pb ■ (2.26) 



The main advantage of the Hamiltonian ( p. 26 ) density with respect to fl2.1|) is that it 
depends on the momenta 7r lJ quadratically and hence the Legendre transformation from 
the Hamiltonian to the Lagrangian description can be easily performed. Explicitly, the 
canonical equations of motion for A, B and g^j take the form 



d t A = {A,H} = v A , d t B = {B,H} =v B , 
2BNk 2 

dtQij = {gij,H} = —g ijkl7r kl + ViNj + VjNi 



It is natural to define 



V9 



K lJ = ^(d t g l ,-V l N J -V J N l ) 



(2.27) 
(2.28) 



so that the equation on the second line in Q2.27 ) implies 

Bk 2 



V9 



(2.29) 



Then it is easy to determine the corresponding Lagrangian 
L = f d^dtg^wij + d t AP A + d t BP B -H) = 



d xN V9 ( -^Ki&WKu - K 2 B(E^g ijkl E kl -A)- k f(A) 



(2.30) 



Finally we eliminate non-dynamical fields A and B from (2.30). Varying ( |2.30| ) with respect 
to B 



~ BW Ki i 0K K ™ ~ ( EiJ Gm E -A)=0 

and varying ( p. 30 ) with respect to A we find 

B-f'(A) = 0. 
Inserting (2.32) into ( 2.31| ) we obtain the equation 



(2.31) 



(2.32) 



/«(A)«a 



Kij&WKu - K 2 (E^G ijkl E kl - A) = 



(2.33) 



-7- 



that, for known function /, allows us to express A as a function of Kij,gij 

A = ^{K ij ,g ij ) . (2.34) 

Inserting these results into the Lagrangian ( |2.30| ) we find the new class of F(R) theories of 
gravity in the form 



L 
C 



d D x£ 



f'{*{Kiji9ij))K 



(2.35) 



Let us consider the specific example of the theory defined by the function f(A) = \/l + A 
1. For this function the equation ( 2.32j ) implies 

1 



B 



2VTTA ' 



Inserting this result into ( 2.31| ) we find 
A = l I' {.(/,, . A', ; . 



ki 



1 - 



(2.36) 



(2.37) 



and hence the Lagrangian density takes the form 



-k 2 N^ a/1 - —K i3 G^K kl J\ + E^G ljkl E kl - 1 



We see that ( [2.38 ) coincides with the Lagrangian density studied in [28]. 



(2.38) 



Finally we show that the action S = f dtL where L is given in fl2.38| ) is invariant 
under the foliation-preserving diffeomorphism that by definition consists a space-time de- 
pendent spatial diffeomorphisms as well as time-dependent time reparameterization. These 
symmetries are generated by infinitesimal transformations 



6x i = x H - x i = C*(t,x) , 5t = t' -t = f(t) . 

It was shown in [12] that the metric transform under ( 2.3S| ) as 

9ij(t',x') = 9ij(t,x) - g u (t,x)djC l - diQ k g kj (t, x) , 
/^',x') = g^{t^) + d n Cg n \t^) + g m (t^)d n ^ 

and the fields iVj(i, x), N{t) transform under d2.39|) as 

NUt',*!) = Ni(t,x) - Ni(t,y)f - N 3 {t,x)diCi - //,,(/. x;C , 
N H (t',x') = iV i (t,x)+^(t,x)5 i C(t,x)-^(t,x)/-C(t,x) 
N'(t') = -N(t)f . 



(2.39) 



(2.40) 



(2.41) 
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Further, the transformation property of the metric components ( |2,40| ) imply following 
transformation prescription for Qijki 

g' ijkl (t',x.') = Gijkiit,^) + diC m {t,x)g mjkl {t,x.) + ^c m (t,x)e? imfc; (i,x) + 
+ g i]ml {t^)d k c m {t^) + g ij km(t,x)di( m (t,x) . 

(2.42) 



It can be also shown that transform covariantly under ( |2.39| ) 

KUt',x.') = Kij(t,x) - K ife (t,x)5,C fc (t,x) - diC k (t^)K kj {t^) 



(2.43) 



Using these formulas it is easy to see that the expression KijQ^ kl K k \ is invariant under 
fl2~39D . 



Finally we consider the expression E l3 Qij k \E . As was shown in [28] E l ' J transforms 



under ( 2.39 ) as 

£^'(t',x') = &(t,x) +d k C(t,y:)E k Ht,x.) + E* h (t,x)d k (i (t,x) . (2.44) 
Then E % ' J 'Gij k iE kl is clearly invariant under ( [2.39D . Finally if we use the fact that 



N'(t')^g(x.',t')dt'd D x.' = N(t)^g(x,t)dtd D x (2.45) 
under ( 2.39j ) we find that the action ( 2.30| ) is invariant under foliation preserving diffeo- 



morphism ( |2.39 ). 



3. The Second Version of f(R) Hofava-Lifshitz Theory of Gravity 

In this section we present the second way how to define f(R) Hofava-Lifshitz theories of 
gravity. We begin with known f{R) theory of gravity with the action 

S = -^J dtdxNy/gf(R^) , (3.1) 

where is four dimensional covariant Ricci scalar. Then in order to find f(R) Hofava- 
Lifshitz theory of gravity we simply replace R^ in ( |3.1[ ) with the combination K^K 1 ^ — 
\ 2 K - E^g ijkl E kl that appears in the Hofava-Lifshitz theory that obeys the detailed 
balance condition 5 . 

As in previous section we introduce the auxiliary fields A and B and write the action 
Q3.1| ) in the equivalent form 

S = j dtd D ^N(B{K tJ K^ - X 2 K 2 - E^Q l]kl E kl - A) + f(A)) . (3.2) 



5 We would like to stress that this procedure can be easily generalized to theories that break the detailed 
balance condition. Generally, we could replace with KijK 1 ^ — \ 2 K — V where the explicit form of V(g) 
was given in pt^ ]. 
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Clearly the solution of the equation of motion for B reproduces the original action while 
the equation of motion for A implies f'(A) = B. Inserting this result into ( |3.2| ) we find 

S = k 2 J dtd D yL^gN{f'{A){K i:i K^ - \ 2 K 2 - E^Q ijkl E k[ - A) + f{A)) . (3.3) 

Our goal is to find the Hamiltonian corresponding to the action ([T^) . Even if this action is 
invariant under foliation preserving diffeomorphism we can consider more general case and 
presume that the theory does not obey the projectability condition so that N depends on 
x. Then we introduce the set of conjugate momenta 7r* J (x), 7r*(x), 7Tjv(x) and Pa(x), Pb(x) 
where 



ir ij = J^- = K 2 J-gBQ^ kl K kl , 



and where the remaining momenta form the primary constraints of the theory 

7T JV (X)«0, 7T i (x)«0, P^xJfsO, Pb(x)«0. 

Then it is straightforward to find corresponding Hamiltonian 

H = J d D xH(x) = J d D x(iV(x)W (x) + N i (x)H i (x) + 



(3.4) 



(3.5) 



H 



o 



+ ^(x)7r i (x) + t. 7 v(x)vr 7V (x) + v A {x)P A {x) + ^(x)P B (x)) 
' M i .2 r^TDiT^iin ipkl , A \ .2 



^Qijkl" + K z ^g-B(E l Jg ijkl E« l + A) - K z ^f(A) 



(3.6) 



Since the constraints vr J (x) ~ , vtat(x) ~ have to be preserved during the time evolution 
of the systems we find 

a t7 r(x) = {ir N (x),H} = -H (x) « , 
d t 7r\x) = {tt<(x),P} = -H l (x) «0 . 

(3.7) 

On the other hand the requirement of the preservation of the constraints Pa (x) ~ , Pb (x) r 
during the time evolution of the system imply 

dtP A (x) = {P A (*),H} = G A (x) « 

(3.8) 

and 

£> t P B (x) = {P B (x), P} = iVG B (x) « , (3.9) 
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where 



G A = -B + f'(A) , 

G B = 1 - Tr ij g m n kl - K 2 ^(E^g ijkl E kl + A) 
K & \/9 



(3.10) 



Now we calculate the Poisson brackets of these constraints when we again introduce the 
notation = (Pa, Ga, Pb, Gb) , J = 1, 2, 3, 4. Then we have 

{$i(x),$ 2 (y)} = -/"(A)(x)J(x-y) , 

{$ 1 (x),$ 3 (y)} = , 

{$i(x),$ 4 (y)} = « 2 J B(x)Vff(x)5(x-y) , 
{<D 2 (x),<I>3(y)} = -*(x-y) , 
{$ 2 (x),$ 4 (y)} = , 

{$ 3 (x),<I>4(y)} = K 2 ^(E^g i3kl E kl + A)(-K)5(-K-y) . 

(3.11) 

We see that 3>/ is collection of the second class constraints that can be set strongly to zero 
and explicitly solved 6 . First of all, <l?i = $3 = imply that Pa = Pb = 0. Further, from 
$2 = we find 

B = f(A) . (3.12) 
Finally the constraint <J>4 = together with ( |3.12| ) gives 

-kl .2 /-/ rnijn r?kl 



^Gijki^ = k yfg(E %:> QijkiE + A) . (3.13) 
Again, for known / the equation above allows us to find A as a function of Gijkl^ 1 



^P(A)Jg 

f 

and y/gE lj GijkiE kl so that 



A = y(JL n i3g ijm kl^ ^E^Gi jk iE kl ) . (3.14) 

Using these results we finally find the Hamiltonian constraint TCq in the form 

2 



(3.15) 



Acknowledgements: This work was supported by the Czech Ministry of Education under 
Contract No. MSM 0021622409. 



6 It the similar way as in previous section we can easily show that the Dirac bracket {gij(x),7r fei (y)} 
coincides with the Poisson bracket {gij(x.),iv kl (y)}. 
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